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ZruSbs is shown to be entropically stabilized versus decomposition to Zn and ZnSb though the 
effects of configurational disorder and phonon free energy. Single phase stability is predicted for a 
range of compositions and temperatures. Retrograde solubility of Zn is predicted on the two-phase 
boundary region between Zn4Sb3 and Zn. The complex temperature dependent solubility can be 
used to explain the variety of nanoparticle formation observed in the system: formation of ZnSb on 
the Sb rich side, Zn on the far Zn rich side and nano-void formation due to Zn precipitates being 
reabsorbed at lower temperatures. 
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A. Introduction 

The Zn-Sb binary phase system has been of interest 
for many years in the search for efficient and low-cost 
thermoelectric materials. Of primary interest has been 
the Zn 4 Sb3 phase which exhibits a thermoelectric figure 
of merit, zT, in excess of 1 in intermediate temperature 
ranges^. This phase, being composed of environmentally 
benign and relatively earth abundant elements, continues 
to draw active research. 

Zn 4 Sb3 exhibits exceptionally low lattice thermal 
conductivity^, due in part to its high configurational dis- 
order. The room temperature structure of Zn 4 Sb 3 (i?3c, 
Figure fil denoted '£-Zn 4 Sb 3 ' here^ 1 ( or sometimes 'e- 
Zn 4 Sb 3 ' 2 ), contains an anionic framework composed of 
30 Sb, divided between 6 Sb 2 ~ dimers and 18 isolated 
Sb 3 ~. With 39 Zn 2+ a charge balanced composition is 
obtained at Zn 13 Sb 10 . Of these Zn cations, at most, 36 
can fit on Zn framework sites and 3 must be distributed 
among 3 crystallographically distinct interstitial sites^. 
For clarity, we will always refer to the entire phase as 
'Zn 4 Sb 3 ' and write 'Zn^Sbio' only when referring to a 
specific configurational composition within 'Zn 4 Sb 3 '. 

Previous ab initio studies of Zn 4 Sb 3 have included 
the electronic density of states and formation energy 
of various atomic configurations ranging in composition 
from Zn^Sbio to Zni 4 SbicP^. These results confirm the 
expectation from charge counting that Zni 3 Sbio has a 
Fermi level in the band gap, Zn-deficient Zn 12 Sbio nas 
a Fermi level in the valence band and Zn-rich Zni 4 Sbio 
has a Fermi level in the conduction band. All reports of 
formation energjEHS also show agreement that all Zn 4 Sb 3 
configurations have a positive formation enthalpy at OK 
with respect to Zn and ZnSb. Yet Zn 4 Sb 3 is observed be- 
low 700K. Furthermore, at low temperature, a reversible 
transition to a low sym metry , meta-stable, ordered 'a' 
(and a') phase is observed^. 

Synthetically, Zn 4 Sb 3 is usually prepared with excess 
Zn: Zni 3 . 3 Sbio- As Zn has a high vapor pressure even 
at the relatively low synthesis and operation tempera- 
ture, as well as the tendency to readily oxidize^, it is 




FIG. 1: (Color online) Conventional hexagonal unit cell of 
Zn4Sb3. The Sb sub-lattice (orange) is composed of isolated 
Sb and Sb dimers. The Zn 'A' sub-lattice (blue) is 90% 
occupied. The 'B' (green), 'C (brown), and 'D' (pink) Zn 
sub-lattices are 5% occupied. The primitive rhombohedral 
cell is 1/3 the volume of the conventional cell. 



not very surprising that excess Zn helps form the phase. 
Zn 4 Sb 3 is always a heavily doped p-type semiconductor 
with room temperature (Hall) carrier concentration be- 
tween 6 — 9xl0 19 /cm 3 |m According to the simple charge 
counting and ab initio calculations, this would be best 
explained by a slight Zn deficiency, Zn 13 _ ( 5Sb 10 , where 5 
should be between 0.016 and 0.024. If n/nu a u is 2 as 
predicted by SinglP, then S should be between 0.008 and 
0.012. It is perhaps surprising that with all the intersti- 
tial Zn sites, synthesis with excess Zn does not produce 
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n-type material. 

In this work, we expand on previous ab initio results 
to consider a thorough thermodynamic investigation of 
Zn4Sb3. A thermodynamic ensemble approach is neces- 
sary because no single configuration (or even small num- 
ber of configurations) can adequately describe the phase. 
By assembling a Grand Canonical Partition Function we 
are able to quantitatively prove entropic stabilization. 
We also predict a region of single phase stability on a 
temperature versus composition diagram that exhibits 
retrograde solubility of Zn with S > that would always 
produce p-type Zn4Sb3 in intermediate temperatures. 

B. Computational Methods 

The natural thermodynamic function to account for a 
range of possible compositions is the Grand Canonical 
PotentiaP(GCP): 

0(T lM ) = -^nfce-(^^)/^ (1) 

where k is Boltzmann's constant, T is temperature, N 
is the total number of atoms in the system, E s is the 
total energy of state s, N s is a vector of the number of 
atoms in each species in state s (with elements summing 
to N) and /i is a vector containing the chemical potential 
of each species. This greatly simplifies for both Zn and 
ZnSb, which do not exhibit configurational disorder: 

(f>{T,^) = — - fJ-- x (2) 
n 

where £o is the energy per unit cell, n is the number of 
atoms per unit cell and xq is the composition in atomic 
fraction. 

For Zn4Sb3, we consider the GCP within an indepen- 
dent cells approximation. We assume that each primitive 
unit cell of 23 ± 2 atoms is non-interacting with neigh- 
boring cells, in the sense that the defect configuration 
present in one cell does not affect the energies of defect 
configurations in a nearby cell. (All our ab initio cal- 
culations are nevertheless performed on infinite periodic 
systems with suitable k-point sampling.) This approach 
is useful in this system because the unit cell is rather 
large. This assumption is validated by computing the 
energy of supercells with different configurations in each 
primitive unit cell and comparing it to that predicted by 
summing the energies of the constituent primitive cells. 
Under the independent cells approximation: 

<S> ( T , /•*) = — ^ /» ' xo 
n 

_ k JLl n /' 1 + ^ m . e -(A ei -M-A„ i )/ fcB T'j (3) 

where £o is the ground state energy per unit cell, n is 
the number of atoms per unit cell and xq is the ground 



state composition. For each configuration, i, mi is the 
symmetric multiplicity, Ae^ is the change in energy from 
the ground state and Arii is the change in the number 
of atoms from the ground state. 

The energy for each configuration is calculated under 
the generalized gradient approximation (GGA) using the 
projector augmented wave (PAW) method with Perdew- 
Burke-Ernzerhof (PBE) potentials as implemented in 
VASP 4.6, neglecting spin orbit coupling. All unit cell 
parameters and atomic positions were allowed to relax to 
find the lowest energy configuration to within 10~ 4 eV. 
A final static calculation was performed for an accurate 
total energy. 

Defect configurations were systematically generated by 
enumerating defect combinations deviating from the un- 
defected Zn^Sbio structure ('A' sub-lattice fully occu- 
pied and 'B','C and 'D' sub-lattices fully unoccupied). 
Defect clusters were composed of a combination of va- 
cancies on the 'A' sub-lattice, and occupation of an in- 
terstitial site on the 'B', 'C or 'D' sub-lattices. Clusters 
of up to 6 defects were enumerated, excluding some struc- 
tures because the defects were too close (nearest neighbor 
A and C sites both occupied) . After allowing the atomic 
configuration to relax, the resulting atomic positions were 
projected onto the closest unrelaxed configuration and 
any duplicate configurations were excluded in order to 
avoid over-counting states. 

Phonon density of states and vibrational free ener- 
gies were calculated using the 'supercell' method as im- 
plemented in the Alloy Theoretic Automated Toolkit 
(ATAT)p2HHl_ Since the computational resources needed 
to compute phonon modes for all Zn4Sb3 configura- 
tions are prohibitive, representative configurations were 
selected at four compositions between Zn 12 Sbio an d 
ZnisSbio configurations. A second Zni3Sbio configura- 
tion was computed as well to assess the error in our ap- 
proximation. For each of these configurations and the 
end-members Zn and ZnSb, phonon modes were calcu- 
lated at 0%, 1% and 2% strain to account for the effects 
of thermal expansion, under the quasi-harmonic approx- 
imation. 

The vibrational contribution to the free energy is incor- 
porated into the GCP through a nested sum in the parti- 
tion function^. For each distinct configuration, phonon 
occupation accounts for small displacements around the 
local energy minimum, resulting in a temperature depen- 
dent free energy correction. 

GCPs were assembled for Zn, Zn4Sb3 and ZnSb. Phase 
equilibrium is determined by equality of two respective 
GCPs. The equilibrium composition of each phase can 
be determined by: 

V fl <j>(T, f i) = -x (4) 

C. Results and Discussion 

Over 100 unique, stable configurations were enumer- 
ated in the /3— Zn4Sb3 primitive rhombohedral unit cell 
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FIG. 2: (Color online) Formation enthalpy at OK for Zn4Sb3 
configurations with respect to Zn and ZnSb. The solid black 
line represents the convex hull of structures confined to the 
Zn4Sb3 lattice. Supercell calculations in green were done us- 
ing the conventional unit cell, each composed of three prim- 
itive unit cells. The dashed line shows the predicted energy 
for supercells following the independent cells approximation. 

with between 21 and 25 atoms per cell (the number of 
Sb was held constant at 10). The OK formation enthalpy 
of these configurations with respect to Zn and ZnSb is 
shown in Figure[2j 

Consistent with previous theoretical studies^^, all con- 
figurations have a positive formation enthalpy with re- 
spect to Zn and ZnSb. The convex hull is plotted in 
solid black and connects the ground state configurations 
(for Zn 4 Sb3 phase) at each composition. The composi- 
tion Zni4Sbio does not have a configuration touching the 
convex hull. This means that for an ensemble of atoms 
of composition Zn 14 Sbio constrained to remain in the 
Zn 4 Sb3 lattice, it would be more energetically favorable 
to form a mixture of cells of composition Zni3Sbio and 
ZnisSbio. If the system is allowed to adopt any lattice, 
then, at OK, it would be even more energetically favor- 
able to separate into Zn and ZnSb crystals. Figure [2] also 
allows us to test the independent cells approximation. 
The data points marked supercells are all made up of 3 
primitive cells of different compositions to yield the in- 
termediate compositions shown. Between Zn^Sbio and 
Zni3Sbio, we see that the supercells produce a slightly 
lower energy by about 3 meV/atom than predicted by the 
independent cells approximation. Between Zn 13 Sbio an d 
Zni4Sbio, the energies from the supercells are either very 
close or slightly higher than predicted by the independent 
cells approximation (the dashed line). The small energy 
difference of a few meV/atom justifies the use of the in- 
dependent cells approximation, although a consideration 
of a possible systematic error is discussed later. 

The formation free energy due to phonons (relative to 
Zn and ZnSb) is shown in Figure[3] for the 5 representa- 
tive configurations. In all cases, the phonon contribution 
to the free energy favors the Zn 4 Sb3 phase over Zn and 
ZnSb. This is likely a result of the softer phonon modes in 
the more complex and open crystal structure. Increasing 
Zn concentration and thereby increasing the Zn-disorder 
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FIG. 3: Phonon contribution to the formation free energy 
with respect to Zn and ZnSb. The two ZnisSbio configu- 
rations are in good agreement to lmeV/atom. There is a 
more favorable contribution to the formation energy for more 
disordered Zn-rich configurations. 



results in a more favorable contribution to the free energy. 
There is good agreement between the two configurations 
of composition ZnisSbio, differing by lmeV/atom at 
1000K. It then seems reasonable to assume that the vi- 
brational free energy of the representative structures may 
be applied to all the configurations at that composition 
when we compute the GCP of Zn 4 Sb3. 

The computed region of single phase stability for 
Zn 4 Sb3 is shown in Figure[4j Including configurational 
and vibrational effects to the free energy, Zn 4 Sb3 is 
found to stabilize at around 700K at a composition of 
Zni2.992SbiQ. As the temperature increases, the range of 
stable single phase compositions increases, more broadly 
on the Zn deficient side. (Eventually, the solid melts, 
but our analysis focuses on the solid-state portion of the 
phase diagram.) Interestingly, on the Zn-rich side, we 
observe retrograde solubility of Zn. That is, the high- 
est stable concentration of Zn in Zn 4 Sb3 decreases with 
increasing temperature. This unusual finding is a result 
of the composition-dependent relationship between the 
enthalpy and entropy of mixing. Typically defects of all 
kinds entropically stabilize stoichiometric variations on 
both sides of a valence precise composition. In the case 
of Zn 4 Sb3, the valence precise composition is already 'de- 
fected' with interstitial Zn. A miscibility gap arises for 
compositions above ZnisSbio due to the large jump in 
formation enthalpy compared to the Zn-dcficient compo- 
sitions. The effect is lessened by the counteractive effects 
of higher vibrational entropy for the Zn-rich composi- 
tions. If the vibrational effects were neglected, the Zn- 
rich boundary would not curve to the right until a much 
higher temperature. 

The region of single phase stability shown in Figure|4] 
has readily observable consequences. The black vertical 
line at composition ZnisSbio represents the valence pre- 
cise structure with a Fermi level inside the band gap. 
At any Zn-deficient composition from Zn 13 Sb 10 , we ex- 
pect a partially filled valence band resulting in a p-type 
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FIG. 4: (Color online) Single phase stability region for 
Zn4Sb3. Solid black line corresponds to the charge balanced 
composition. Stable compositions to the left of the line will 
result in a p-type semiconductor. Inset shows a close-up 
of the region of retrograde Zn solubility. The orange band 
corresponds to compositions estimated from charge carrier 
measurements^. 

semiconductor. (Surprisingly, there is a recent report of a 
Zni 3 Sbio configuration being slightly p-type^. Our elec- 
tronic density of states results are in agreement with pre- 
vious studies^].) All stable compositions below 1000K 
in Figure[4] result in a p-type semi-conductor. Attempt- 
ing to dope Zn4Sb3 with excess Zn will not result in an 
n-type semiconductor but a two-phase equilibrium be- 
tween Zn metal and p-type Zii4Sb3. Furthermore, with 
the consideration of retrograde solubility, if Zn-saturated 
single phase Zn 4 Sb3 is heated up, it passes through a re- 
gion where it becomes thermodynamically more stable to 
precipitate Zn metal. With more heating, the Zn is re- 
absorbed as more Zn-rich compositions are stable. Upon 
cooling, the same precipitation and absorption should oc- 
cur if held in thermodynamic equilibrium. 

On the Zn-deficient side of the stable phase region the 
usual temperature dependent solubility is observed. Here 
one would expect ZnSb to precipitate as the temperature 
is reduced. Because this would occur at low tempera- 
ture in the solid state, small nanometer sized precipitates 
would be expected. Such particles are indeed observed 
in Sb-rich samples^. 

The high temperature phase boundary on the Zn-rich 
side could explain the formation of Zn nanoparticles ob- 
served in some Zn4Sb3 samples^. Cooling a Zn-rich com- 
position above Zn 13 Sbio would precipitate Zn (possibly 
nano-particles if cooled fast enough 19 ). Upon further 
cooling through the retrograde region, some of the Zn 
would be reabsorbed into the /?— phase. The absorption 
of nanoparticle Zn may explain the nano-voids observed 
in some sampled^] 

Experimental carrier concentration measurements 
from Toberei^ from single phase (3— Zn4Sb3 samples 
range from 6 — 9xl0 19 /cm 3 . This corresponds to a com- 
position range of 0.002% atomic Zn between Zn^.gssSbio 
and Zn 12 .992Sb 1 o assuming n/nuaii = 2 P. This concen- 
tration range, corresponding to the orange band in the 



inset of Figure[4j is in close agreement with the predicted 
stable compositions near the stabilization temperature. 
This composition range is much smaller than the 0.2% 
atomic Zn observed by microprobe analysis. This dis- 
crepancy is not surprising since the composition range is 
basically the limit of the microprobe resolution. 

There is some expected uncertainty in our predicted 
temperature of stabilization of Zn4Sb3. Our predicted 
temperature of 700K is much higher than might be ex- 
pected. There are several possible explanations for this. 
One could be a deficiency in our methodology to fully 
account for the sources of entropy in such a complex dis- 
ordered structure. Our enumeration method may under- 
estimate the number of Zn4Sb3 configurations by com- 
paring the relaxed structures to the fixed lattice in Fig- 
ure[lj which represents, in a sense, a measured average 
over many possible configurations. The vibrational con- 
tribution to the free energy has a significant effect in 
lowering the stabilization temperature. Assuming a rep- 
resentative structure for each composition might be an 
underestimate for the more disordered configurations and 
the effect could be even more pronounced. Lastly, we 
noted earlier that the supercells on the Zn-deficient side 
of Zni3Sbio had slightly lower formation energy than ex- 
pected under the independent cells approximation. This 
neglected cell-to-cell interaction ( 2meV/atom) could 
lower the stabilization temperature by several hundred 
degrees. In either case, we expect the retrograde solubil- 
ity to remain, and possibly intensify, because the super- 
cells on the Zn-rich side of Zni3Sbio are in good agree- 
ment with the independent cells approximation. 

To ensure that our finding of a single-phase region for 
the /3-Zn4Sb3 lattice is not an artifact of the independent 
cell approximation, we have conducted separate Monte 
Carlo simulations based on a cluster expansion Hamilto- 
nian fitted to our database of structural energies. These 
simulations confirm the presence of a single phase region 
over the temperature range where /3-Zn4Sb 3 is stable. 
Hence, the independent cell approximation was deemed 
reliable and was used throughout this work. It provides 
a convenient explicit expression for the free energy and is 
immune to the fitting errors inherent to the Hamiltonian 
construction procedure. 

Finally, we consider how Zn4Sb3 might interact with 
other nearby phases, namely, a— Zn4Sb3 and ZnsSby. 
First principles calculation of a— Zn 4 Sb 3 yields positive 
formation energy of 19meV/atom with respect to Zn 
and ZnSb. This is 3meV/atom below the lowest energy 
ft— Zn4Sb3 configuration of composition Zni3Sbio- As- 
suming no configurational disorder in a— Zn 4 Sb3 we pre- 
dict j3— Zn4Sb 3 to become energetically favorable com- 
parted to a— Zn4Sb 3 at 300K (although both are still 
meta-stable with respect to Zn and ZnSb at this tem- 
perature). It seems possible that there is a temperature 
range in which meta-stable P—Zn^Sh^ is observed before 
becoming thermodynamically stable with respect to Zn 
and ZnSb at some higher temperature. 

Another phase of recent interest is ZngSby, which has 
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FIG. 5: (Color online) Possible phase stability regions of 
Zn4Sb3 and ZngSby under conditions favorable to the forma- 
tion of ZngSby. The dotted line represents the phase bound- 
ary of Zn4Sb3 in the absence of ZngSby. 



been characterized experimentally and studied through 
first principles calculations^^!!. These calculations re- 
veal that, even though ZnsSby is unstable in bulk form, 
it could be stabilized in nanocrystaline form if surface 
stress or energy contributions result in a free energy de- 
crease of 5mcV/atom (relative to ZnSb and Zn 4 Sb3). 
Taking these assumptions into effect, the resulting phase 
diagram is shown in Figure|5j We see that Zn 8 Sb7 would 
stabilize at a temperature higher than Zn 4 Sb3 and a two 
phase region results. This also reduces the stable Zn- 
deficient Zn 4 Sb3 composition range predicted in the ab- 



sence of Zn§Sb7 (dashed line in Figure[5]). Other uniden- 
tified phases as reported in some phase diagrams^ could 
have a similar effect. 



D. Conclusions 

From our first principles investigation, we have shown 
that Zn 4 Sb3 is entropically stabilized with the help of 
configurational and vibrational entropy. Under the in- 
dependent cells approximation of the grand canonical 
potential we predict a region of single phase stability 
near Zn 12 .992Sb 10 , which results in a nominally p-type 
semi-conductor. Additionally, we predict a temperature 
range with retrograde Zn solubility. The temperature de- 
pendent solubility can be used to explain the variety of 
nanoparticle formation observed in the system: forma- 
tion of ZnSb on the Sb-rich side, Zn on the far Zn-rich 
side and nano-void formation due to Zn precipitates be- 
ing reabsorbed at lower temperatures. 
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